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HEAT TRANSFER IN ISOTROPIC TURBULENCE 
DURING THE FINAL PERIOD OF DECAY 
By D. ¥. Dunn and ¥. H. Reid 

SUMMARY 


The problem of heat transfer in isotropic turbulence with a constant 
mean temperature gradient is considered during the final period of decay. 
The Reynolds and Pdclet numbers are then very small, and all triple cor- 
relation terms can be neglected in the equations for the double correla- 
tions. On this basis, it is found that the temperature field ultimately 
becomes independent of the initial conditions on the temperature and has 
characteristics determined only by the mean temperature gradient, the 
physical properties of the fluid, and the characteristics of the turbu- 
lence. Detailed analytical and numerical results are obtained for the 
asymptotic state. 

The mean turbulent heat transfer is in the direction of the mean 
temperature gradient, with a magnitude proportional to the magnitude of 
the latter. Although it approaches zero when the Prandtl number approaches 
zero, its dependence on the Prandtl number Is not large for Prandtl num- 
bers of order unity and larger. This type of Prandtl number dependence 
is typical for many of the other results depending on both velocity and 
temperature fluctuations. In fact, the rate of decrease with separation 
distance of the two-point temperature-velocity correlation varies little 
over the full range of Prandtl numbers and is always about the same as it 
Is for the double velocity correlation. In contrast, all results involving 
only temperature fluctuations display a strong dependence on the Prandtl 
number. For example, for small Prandtl numbers the double temperature 
correlation falls off much more slowly with separation distance than the 
velocity correlation does, while for large Prandtl numbers the opposite 
Is true. 


INTRODUCTION 


The simplest case of turbulent heat transfer is the problem first 
considered by Corrsin (ref. l) , in which the temperature of the fluid is 
specified to have a constant mean gradient in some preferred direction, 
while the velocity field, which is assumed independent of the temperature 
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field, is regarded as isotropic and known. Thus, though mean values 
associated with the velocity field only are isotropic, those associated 
with the temperature field are axisymmetric . A study of this problem, 
although it is highly idealized, is expected to give some idea of the 
nature of turbulent heat transfer and, by analogy, also of turbulent mass 
transfer . 

Such a study is also a natural preliminary step before attempting 
the much more complicated problem of turbulent momentum transfer. The 
simplest example of the latter, the so-called homogeneous shear flow with 
a constant mean velocity gradient studied by Reiss (ref. 2) and later by 
Burgers and Mitchner (ref. 3)> is mathematically similar to the present 
problem in some respects, though with additional complications due to 
the vector nature of the transported quantity. The present problem, 
because of the scalar nature of the transported quantity and the resulting 
more restrictive symmetry conditions, is considerably simpler mathemati- 
cally. Its solution may be expected to provide useful clues towards the 
solution of the more difficult problem of homogeneous shear flow. 

In addition to being of theoretical interest, this idealized problem 
of turbulent heat transfer is of considerable practical interest, since 
a direct experimental investigation of it appears to be feasible. For 
example, if the horizontal bars of the usual turbulence-producing grid 
in a wind tunnel are heated with a proper variation in the temperatures 
of successive bars, the conditions of the idealized problem should be 
reproduced fairly closely in a region downstream of the grid. There are, 
of course, various approximations involved in relating the experimental 
situation to the theoretical problem, but they appear to be no more serious 
than the ones normally associated with the wind-tunnel study of Isotropic 
turbulence . 

Although the problem as stated appears to be the simplest problem 
of turbulent transfer that can be formulated, there are still formidable 
mathematical difficulties present in the general case, when the inertial 
transfer terms are dominant in the basic equations. In particular, the 
system of equations for temperature mean values that can be derived from 
the heat-transfer equation is incomplete, just as is the well-known system 
for velocity mean values derived from the dynamical equations. Of course, 
as in the latter case, one could, to obtain a determinate system, mahe 
some type of assumption such as the one most commonly made at the present 
time that fourth-order mean values are related to second-order ones as 
they would be for a normal probability distribution. 

Froudman and Reid (ref. 4) have recently used this statistical 
assumption in an analysis of isotropic turbulence involving mean values 
of velocities at several pointB but at the same time (see also ref. 5 
by Tatsumi) . Reid (ref. 6) has applied the idea in a similar manner to 
the study of the convective effects of isotropic turbulence on arbitrary 
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vector and. scalar fields, including, in particular, the case of isotropic 
temperature fluctuations in isotropic turbulence. Chandrasekhar (ref. 7 ) 
has used it in an alternative approach to isotropic turbulence, in which 
mean values at different times as well as at different points are con- 
sidered. At the present time, however, the possibilities of the assump- 
tion have not been fully exploited even in these problems involving com- 
plete isotropy. Until further progress has been made with the latter, it 
does not appear worthwhile as yet to pursue this approach further in the 
present more complicated problem. 

The objectives of the present study, as well as of the earlier inves- 
tigations of Corrsin (ref. 1) and Mdlfese (ref. 8), are considerably more 
limited, being restricted to results that can be derived without a detailed 
consideration of the triple and higher order correlations. Corrsin and 
Mdlfese have obtained in this manner results of a general nature for a 
variety of different situations but have not treated any special case in’ 
detail. In the present investigation, attention will be restricted mainly 
to the final period of decay, and for this case essentially complete 
results will be derived for the double temperature correlation and the 
double temperature-velocity correlation and their Fourier transforms. 

During the final period, the turbulence Reynolds numbers are so low that 
inertial transfer effects are very small, and the triple-correlation terms 
in the equations for the double correlations can be neglected. To this 
order of approximation the problem is determinate, and the statistical 
hypothesis mentioned above is not required. The conclusions may be 
expected to give an approximate description of the decay of the velocity 
and temperature fluctuations far downstream of the grid in the experiment 

1 described earlier. 

, Although some of the most interesting features of turbulent heat and 

momentum transport are obviously absent during the final period of decay, 
the latter must still be recognized as an important aspect of the general 
problem. For this reason, and because of the relative simplicity of the 
analysis, a detailed treatment of this special case is appropriate first, 
before proceeding with the more difficult general case. 

This research was carried out at the Department of Aeronautics, The 
Johns Hopkins University, under the sponsorship and with the financial 
assistance of the National Advisory Committee for Aeronautics, and at the 
Ballistic Research Laboratories, U. S. Army, Aberdeen Proving Ground. The 
authors would like to thank Prof. Stanley Corrsin for many helpful dis- 
cussions. Blanks are also due to Mrs. Frances Chenoweth and Mrs. Barbara 
Jacobs for assistance with the numerical calculations. 
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SYMBOLS 

r 

E(k) 

defining scalar of isotropic energy spectrum 

function 

f 

f(r) 

defining scalar of Rij^rj; longitudinal velocity cor- 
relation coefficient 

... 

G 

tK» 

vector with components, G^ = 

cix^ 

. — 

k 

wave -number vector 

- 

L(r,m) 

defining scalar of Lj^r,^ 


L'(r,m) 

_ , ih 

defining scalar of L-j/^r,Aj 
7—7 7 , . . .1 



= §[^(£> t ) u i(x + ^ t ) + ^(x + jC^t)\i i (x,t)J (axisymmetric) 


L l'(z>h) = 5|^fe t ) u i(x + £jt) - d(£ + frtju^j^t)] (axisymmetric) 



thermal integral scale (see eq. (87)) 

4 

M(r,m) 

temperature correlation coefficient (see eq. (96)) 

H 

m. (Ml 

r 

= cos a 


N(r,m) 

temperature -velocity correlation coefficient (see 

eq. (Ill)) 


% 

P£clet number, (u?) ^ A $/k, 


P 

pressure 


% 

heat-transfer correlation coefficient (see eq. (77)) 


R ij(£) = ^ 

L i(x,t^Uj^x + r,t) (isotropic) 

■ 

r a 

Z— 4/2 

turbulence Reynolds number, \u^ J A/v 

0 



NACA IN 4186 


5 


Z 

T 

t 

U P 

^2 


separation vector in correlation functions 

temperature 

time 

initial time 

component of u in a direction normal to A 

<v 

component of u parallel to A, . UjAj_ 
velocity vector with components u^fx^t) 


u 

/v 




x i 

a 

P 


7 



space coordinate 

angle between vector r and unit vector A 

magnitude of mean temperature gradient (a positive 
constant), |vp| 

angle between vector k and unit vector A 

A/ /V 

Kronecker delta 
permutation symbol 


0 (fe'£) 3 Q(k, u) 

2itk 2 atk 2 
0(k) 


Fourier transform of 0^r,Ay (axisymmetric) 
three-dimensional temperature spectrum function. 



9(r,A) = 6(r,m) 3(x,t)-0(x + r,t) (axisymmetric) 
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K 

A(k,|i) 

A'(k,|i) 

A i(k A ) 

V(*,a) 

A(k) 


deviation of temperature from its local mean value 
thermal diffusivity 
defining scalar of A^^k,A) 

defining scalar of A^'^A) 

Fourier transform of Lj_(r,^) (axisymmetric) 

Fourier transform of Li'(r,A) (axisymmetric) 

three-dimensional heat-transfer spectrum function. 



A 


Ai = 


\vi | 


turbulence microscale (see eg.. (67)) 


A 


thermal microscale (see eq. (82)) 

unit vector in direction of mean temperature gradient 
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v 

P 

a 

*u(s) 


0) 


kinematic viscosity 
density 

Prandtl number, v/ic 

Fourier transform of Rj.ifr'} (isotropic) 


vorticity vector with components 


= 


^2 
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Subscripts : 

i,^k ,1 

m 

n 

o 

P 


tensor indices 

maximum 

normal 

at time t Q 

parallel 


FORMULATION OF PROBLEM 


For the problem under consideration, it is assumed that the velocity- 
field is isotropic and unaffected by temperature variations. The gov- 
erning equations for the velocity and pressure are then the Navier-Stokes 
equation 


dui 

dt 


+ 


u^ 


8uj 

8xj 


- — + vV'TIi 

p 1 


(1) 


where V 2 Is the Laplacian operator 




, and the continuity- equation 


du-j 

8x^ 


= 0 


( 2 ) 


The quantities p and v are the density and kinematic viscosity, respec- 
tively, and are assumed to be constant. ,It is supposed that in the fluid 
a constant mean temperature gradient (3 Is maintained in the direction of 
a unit vector A by some external agency (± . e . , £ = | tfF | and 


Ai = 


dT/dx^' 

m j 


so that 


T(x,t) = T c + 3A iXi + d(x,t) (3) 

where T 0 is a constant and d^x,t) denotes the deviation of the tem- 
perature from its local mean value. From the heat-transfer equation 
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the governing equation for the temperature fluctuations is obtained as 

I! + ^7 + ^ (5) 

3 

where k is the thermal diffusivity, also assumed to be constant. The 
equation satisfied by the mean temperature, which is obtained by averaging 

8iSUa 

equation (4), reduces in this problem to the obvious result that - = 0 

dxj 

that is, that the mean turbulent heat transfer vector is a constant. The 
temperature field, because of the preferred direction A, is seen to be 

axisymmetric rather than isotropic. 

The preceding formulation evidently also applies to turbulent mass 
transfer, if the temperature T is replaced by concentration (assumed 
small), and the thermal diffusivity k is replaced by the diffusion 
coefficient. 


The various physical restrictions, such as restrictions on the 
magnitudes of the velocity and temperature variations, which are required 
in order for the above formulation to give an adequate approximation to 
reality, are well-known (see Corrsin's comments in refs. 1 and 9, e.g.). 


During the final period of decay when the turbulence Reynolds num- 
ber is small, molecular transfer effects are expected to dominate inertial 
transfer effects. Quantitative estimates of their relative importance 
can be obtained by order-of -magnitude considerations. For example , in 
the dynamical equation (l) the viscous terms vV^u^ are in this case 

among the dominant terms, and the inertial transfer terms Uj dujJ&xj 


and pressure term -(l/p) (dp/Bx^) are of relative order of magnitude 
/ — n 1/2 /_sl/2 

Vu 2 J A/v, where A is the microscale. The quantities Vu^/ 


and A 


evidently represent a significant velocity and a significant length, 
respectively, for the energy-containing eddies of the turbulence as well 
as for the dissipative eddies under these circumstances (i.e., the micro- 
scale is now of the same order of magnitude as the integral scale) . Sim- 
ilarly, in the heat-transfer equation ( 5 ) the thermal- conduct ion terms 


are among the dominant terms, and the inertial transfer terms 
dxj are of relative order of magnitude q(u 2 Jv, where a 

is the Prandtl number and A 3 is some length characteristic of the 



temperature fluctuations. The thermal microscale, defined by Corrsin 
(ref. 9 ) for the temperature field in a manner analogous to the defini- 
tion of A for the velocity field, is evidently a significant choice 
for the length A 3 . 



NACA TN if 186 


9 


Hence, the final period of decay is characterized by the inertial 
transfer terms in equations (l) and (5) being negligible or more precisely 
by the satisfaction of the following conditions: 


(^) 1/2 x 


« 1 


,-, 1/2 

V 2 ) 

— — « 1 


The first condition implies that the Reynolds number must be small, as 
was to be expected, and the second implies that the Pdclet n umb er must 
also be small. When the Rrandtl number is not too large, it is reason- 
able to expect that A^ and A are of the same order of magnitude, so 

that the Pdclet number is of the order of the Reynolds number and the 
two conditions are essentially equivalent. 

In all subsequent considerations, the approximations described 
above will be accepted, and the dynamical equation and heat -transfer 
equation will be used in the approximate forms 


dt 1 


|i + p*juj = «£* 


The relative errors in these equations are O^R^J In the first and 
O^P^j in the second, where R^ and P^ are the Reynolds number 
and Pdclet number defined in relations ( 6 ) and (7 ) , respectively. 

The situation under consideration, then, is as follows. At some 
time t Q after the initial instant of generation of the velocity and 
temperature fields, the turbulence intensity has decreased to such an 
extent that the Reynolds and Pdclet numbers are both very small and 
equations ( 8 ) and (9) are satisfactory approximations. With their use, 
the subsequent variation of the velocity and temperature fluctuations 
can then be determined to the same order of approximation in terms of 
conditions at time t 0 . 
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In general, of course, the relations between conditions at time to 
and those at the initial instant of generation of the velocity-temperature 
field are beyond the scope of this approximate analysis. These relations 
depend in a complicated and unknown manner on the neglected terms in equa- 
tions (8) and ( 9 ), which are initially important if the initial turbulence 
level is high enough. For present purposes, the conditions at time t Q 

are assumed to be known, at least in a statistical sense. The considera- 
tions will be limited mainly to the double temperature correlation dV 
and the temperature -velocity correlation ■6Uj_'. The results to be pre- 
sented for the double velocity correlation u^Uj ' are well known, but 
they will be included for the sake of completeness. 


Solutions of Equations for Fourier Transforms of 
Correlation Functions 

Let the correlation functions be defined as 



+ i. 



1 

2 


5 (*) u i(* + £) 


${ x + r 


Ms) 


Ll '(~'~) “ 2 S (i) u i(£ + ~) ’ * ~)M~) 


( 10 ) 


(IX) 

( 12 ) 


where the dependence on A has been indicated explicitly to emphasize 

the axisymmetric character of these mean values. For the sake of brevity, 
the dependence on the time will not be indicated. It is to be noted that 
mean values as well as fluctuating quantities are time dependent, since 
a decaying velocity-temperature field is being considered. With the aid 
of equations (8) and (9), the correlation functions, equations (10) 
to (12), are seen to satisfy the equations 


+ 2pA i L i = 2 kV 2 0 
ot 


( 13 ) 
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8Lj 

8t 


+ pAjRij = (k + v)^ 2 !! 


(14) 


= (* + v)V 2 L t 


(15) 


where = u-j^x^uj^x + rj is the double-velocity-correlation tensor 

us ually denoted by this symbol, which to the present order of approxima- 
tion satisfies the equation 


— ^ = 2vA7% i 1 
St 1J 


(16) 


Equations (13) to (15) are special cases of the more general equations 
derived by Chandrasekhar (ref- 10) . 

Since the velocity field is assumed to be isotropic, Rij(r) is an 

isotropic tensor, solenoidal in both of its indices. It can therefore 
be represented in the usual way in terms of the longitudinal velocity- 

correlation function u^f(r), where u2 denotes the mean-square value 
of one velocity component. On the other hand, the vector Li^r,A) is 

axisymmetric, so that in general it requires two defining scalars and 
can be expressed in the form 


L i(£’&) = L l r i + L 2*i 


(17) 


where Lj_ and Lp 


are arbitrary functions of 


r.r = r2 


> 


(18) 


r.A = rm 


) 
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When Li^r^7\j is solenoidal, as it is here, it requires only one defining 

scalar L = L(r,m), which, as Chandrasekhar (ref. 11) has shown, is related 
to L]_ and L 2 in the manner 


1*1 _ -^rmDj. + 

(19) 

L 2 = ^r 2 !^ + rmDm + 2^L 

(20) 


where the differential operators D r 


and Djq have the meaning 



_S_ 

dr 



_d_ 

dm 


m 2 d 
r 2 dm 

> 


( 21 ) 


Exactly analogous considerations apply to the vector Lj_' and its 
defining scalars Iq/ and L 2 ' , which can be expressed in terms of a 
single scalar L 1 as described by relations (19) and (20). 

Equations ( 13 ) to ( 16 ) could be reduced to equations for the defining 
scalars described above, if desired. This is the procedure followed by 
Chandrasekhar (ref. 10). In practice, however, the desired information 
can most easily be derived by considering the Fourier transforms of these 
equations . 

Thus, the Fourier transform of 6^r, is defined in the form 


Q(k,n) 

2jtk 2 


(2 *)' 5 / e (£^) e ' i ~‘~ 


dr 


( 22 ) 


where k*A = kp and where 


©(k,p) has been so normalized that 
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p“ pi 

= / / e(k,n) dp dk (23) 

JqJ -2. 


Likewise, the Fourier transforms of Li^r,A 
defined in the manner 


and Lj_'^r,A^ can be 


A ± (k,A) = (2rt)" 5 dr (24) 

with a similar relation for the transform A^' of L^'. Finally, equa- 
tions ( 13 ) to (l6) can be transformed to 


and 


^ + 4 itpk 2 A^Ai = -2kk% 


(25) 


8Aj_ 

8t 


+ pAj'&.y = -(« + v)k 2 Ai 


(26) 


9Ai' o 

— -(k + v)k ii A i ' 


(27) 


8$ 


8t 


1,3 = -avk^ij 


( 28 ) 


where is the Fourier transform of in the sense of 

equation (24). 

Since $ij(^k) is an isotropic tensor it can be written in the well- 
known form 
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®u(£) 



k l k j\E(k) 
k 2 J 4itk 2 


(29) 


vhere E(k) is the isotropic energy spectrum function. The representa- 

( Sid \ 

tion of the orthogonal fi.e., = 0, corresponding to Oj axi- 

V ar i / 

symmetric vector A^^k,^ in terms of a single defining scalar can, by 

the technique of Proudman and Reid (ref. 4 ), be written directly in the 
form 




A(k,p.)Aa 


(30) 


where A is a function of k-k = k 2 and k-A = kp. Similarly, there 

is an analogous expression for the orthogonal axisymmetric vector A-j_'. 
in terms of a single scalar function A' . Thus, the scalar forms of 
equations (25) to (28) become 


a© 

at 


+ 


2*k 2 0 = -4itp (l 


H 2 )k 2 A 


(31) 


^ + (* + v)k 2 A = -p 

at 


E 

4 rtk S 


(32) 


aA 1 

dt 


+ 


(k + v)k 2 A' = 0 


(33) 


— + 2vk2E .= 0 (54) 

at 
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The derivation of the solutions of these equations satisfying given 
Initial conditions at time t© is straightforward, and the results are 
as follow: 


E(k,t) = E 0 (k)e 


-2vk 2 (t-to) 


( 35 ) 


A(k, p,t) = - 


pE 0 (k) 


4it(K - v)k^ 


-2vk 2 (t-t Q ) _ e -(K+v)k 2 (t-t 0 ) 


Ao(k,ti)e 


-( K +v)k 2 (t-t Q ) 


( 36 ) 


A'Ck^t) = A 0 , (^^)e- (,t+v)Jc2(t - t ° ) 


( 37 ) 


e(k,ti,t) = 


P 2 (l - p 2 )E 0 (k) 


2(k - v) 2 ^ L 


2 


-vk^(t-t 0 ) -Kk 2 (t-t D ) 


- e 


4-jtpfl - p 2 ) 

0 o (k,Ki) + Ao(k,[x) 


K - V 


-2Kk 2 (t-t D ) 


4rtp(l - II 2 ) ^(^^^(K+vJk^t-to) 


K - V 


( 38 ) 


All variables have been explicitly indicated for emphasis, and E 0 (k) = 
E(k,t 0 ), Ao(k,p) = A(k,p,t 0 ), and so forth. Erom the discussion in the 
previous section, it may be expected that for t ^ t Q the errors in these 
results (as compared with the exact results) are i* 1 the first and 

of EH in the others. 
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Asymptotic Behavior of Solutions for Large 
Values of Time 

In the preceding solutions there is still a considerable degree of 
indeterminacy due to the presence of the functions E 0 , 0 O , Ao, and Aq'. 

These express the ultimate influence of the past history of the velocity- 
temperature field during the period of time from the initial instant of 
generation to the time t Q , when the inertial and pressure effects are 
significant. For sufficiently large values of t - t Q , however, it will 

be shown that the results given by equations (55) to ( 58 ) take on much 
simpler forms asymptotically, which depend only to a slight extent on 
the previous history. The asymptotic forms turn out to depend critically 
on the behavior of the functions E, 0, A, and A' for small values of 
k and thus on the behavior of the corresponding correlation functions 
for large values of r. A complete analysis of these questions is com- 
plicated and lengthy and requires the use of the exact equations of the 
velocity and temperature fields. Detailed answers are available only 
for the velocity field alone, in the recent work of Batchelor and Proudman 
(ref. 12) . 

The results of reference 12 indicate that for small values of k 
the behavior of E is in general given by 

E(k,t) = C(t)k^ + o(k5iog e k) (59) 

On the other hand, the exponential factor in equation (55) becomes very 
small as t - to for all except very small values of k, so that the 

asymptotic form is obtained by replacing E Q = E(k,t 0 ) by its leading 
term when k is small. Thus, the asymptotic form for E as t - t Q -* 60 
is 


E(k,t) « C c k 


k -2vk 2 (t-to) 


(40) 


where C 0 = c(t 0 ) . This is the well-known result for the energy spectrum 
function E during the final period of decay. 

In the present problem of the combined temperature-velocity field, 
similar results can be obtained for 0, A, and A’ , providing that plau- 
sible assumptions are introduced at an appropriate stage in the analysis. 
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Equations (22) and (24) are first considered and put into more usable 
forms by the introduction of spherical coordinates (r, a, p) , with 
the aid of the definitions 


r*r = r^ 




r.A = r cos a = rm 


k*A = k cos 7 = kp. 

£p-£p = 


> (U) 


&a-£n 



cos p 




That is, kp and rp are the projections of the vectors k and r 
along A, and &l and rj! are their projections in a plane normal to A, 
so that k = kp + kQ and r = rp + r^. With the aid of the formula 
(ref. 13 , p. 87 ) 



e -iA cos 


P dp = 2 kJ q (a) 


(42) 


where J Q is the Bessel function of the first kind of order zero, the 

integration with respect to the angle p can be immediately carried out, 
leading to the results 


e(k,tx) 

2 *k 2 



cos(krmp.) 




dm dr 



18 


NACA TN 4186 


A(k,p) = — f r r‘^Lp(r,m) cos(krmp) J kr l/(l - m 2 )(l - p 2 ) 

i - v. 2 Jo J-l °L 


dm dr 

m 


A ' (4,p) « — f P r 2 ! ' (r,m) sin(krmp) J h 

1 - \± 2 J 0 J -1 ^ L 


—2 r* CO nl 


“ 

-j 


’ 2 2) 

1 - m jyi - p ) 


dm dr 


( 45 ) 


In these relations, 


Lp — AjLi 


> 


Lp’ = AiLi* 


(46) 


and use has been made of the fact that 0 and Lp are even functions 

of m and Lp' is an odd function of m. These properties are easily 

derived from relations (10) to (12) and (17) • Thus, both 0 and A 
are real and even in p, and A 1 is purely imaginary and odd in p. 

On a purely formal basis for the moment, the first few deriva- 
tives with respect to k of ©, A, and A' at k = 0 are now deter- 
mined from relations (43) to (45). In the case of ©, the results are 
immediately seen to be 
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% 


« 


= 2 ( 2 rt ) _1 r f r^Cr^m) dm dr ( 48 ) 

8k J 0^-l 


For the- quantities A and A', the following relation ( indicated by 
Mdlkse, ref. 8 ) is required: 


L P = 


1 - m2 8 r 3j. 8 

m(l - in 2 ) L 

r 2 8r 8m 

\ I 


(49) 


A similar relation for Lp' follows from equations (17) to (21). With 
the use of equation (49) , and with the assumptions that r^L -» 0 and 
r3L> -» 0 as r -*«>, the following results are obtained: 


A(0,p) = 0 


(50) 


f^(0, p) = 0 

ok 


(5i: 


2 

a A f 


8 k" 


0,\x) = - — r f r i| Lp(r,m) m2p2 + i(l - m 2 ) (l - p 2 ) 

1 _ ,,2w 0 w _i * 2^ ! ' 


1 - p‘ 


dm dr ( 52 ] 


A'(0,p) = 0 


(53) 


8a\ 


8 k 


0,p) = -i H f P r^niLn ' (r,m) dm dr 

1 - p 2 J 0 J -1 


(54 


(Relations ( 50 ) and (53) were obtained previously by Mdlfese.) Hence, the 
behavior of 8 , A, and A', respectively, for small values of k is 
given by 
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e(k,p,t) = c’(t)k 2 + o(k 2 ) 

(55) 

r 

A(k,[j.,t) = C"(p>t)k 2 + o(k 2 ) 

o\ 


A'(k,(j.,t) = C'"(ij.,t)k + o(k) 

(57) 



where 



C"(jx,t) 




(58) 


t 


(See eqs. (48), (52), and (54).) 

The formal results derived above will be valid if the integrals 
involved in equations (47) to (54) actually exist and if the assumptions 

r^L -» 0 and r5L ' -» 0 as r -> <» are valid. For the purposes of the 
present investigation, these conditions will be assumed to be satisfied. 

As mentioned previously, a complete justification would require an elab- 
orate analysis of the behavior of the correlations of temperature and 
velocity for large values of r, after the manner of the work of Batchelor 
and Proudman (ref. 12). Although such an investigation has not been car- 
ried out in full, a partial study reveals no essentially new difficulties 

^Sufficient conditions for the existence of the integrals are that 
r5Lp -» 0 and r^Lp' -> 0 as r -» ». These are not necessary conditions, 
however, and it is assuming less just to assume that the integrals exist. 
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and tends to support the assumptions stated. The latter, moreover, appear 
quite reasonable from an Intuitive point of view. 

By the same type of reasoning as that employed in deriving rela- 
tion (40) , it now follows that the asymptotic forms of relations ( 36 ) 
to ( 38 ) for large values of t - t 0 are obtained by replacing Ep, Ap, 

Aq', and 0 O by their leading terms when k is small, as given by equa- 
tions ( 39 ) and ( 55 ) to (57)* Hence, for t - t Q -» «, relation ( 36 ), for 
example, becomes 

A (k,n,t) » — [ e - 2 vk 2 (t-t 0 ) _ e -(K+v)k2(t-t 0 )l + 

4jt(k - v) u 

C Q " ( (x) k^e “ C v ) k^C t-t Q ) (59) 


where C Q = C(t 0 ) and Cp"(u) = C"(n,tp). Analogous results follow for 
relations (37) and ( 38 )- A further study of these new formulas now shows 
that actually the leading terms are dominant. Consider equation (59), 


for example, and introduce the nondimens ional variable x 
and the reference value 


k\2v(t - t G ) 


An - 


PCp 

4it(« - v) 


/ 2 1+cr „ 2\ 

V 1 - e 1 ) 


( 60 ) 


where a = v/k is the Prandtl number, and xq is defined by 


(20 . 1 + ff\ 

X 1 = I,— l0Se —) 


1/2 


(61) 


That is, Am is the maximum value of the first term in A. Then equa- 
tion (59) can be written as 
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_A 

Am 




- e 


i ±£ X 2 

2a 



+ 


Cq W (h) 

2A m v(t - t Q ) 


x2e 


1+CX o 

- — x^ 

2 a 


( 62 ) 


Since the second term is 
totic form of A when t 


o[l/(t - t 0 )] , the final result for the asymp- 
- to -* 00 is seen to be 


A(k,n,t) « 


PC 0 


4*(k - v) L 


,-2vk 2 (t-t 0 ) 


e -(ic+v)k 2 (t-to) 


(65) 


Similarly, if 0^ is defined as the maximum value of the first term 

in the expression for 0, it is found that the remaining terms in the 
expression for ©/©m are 0(l/(t - t 0 ) ) . Thus, the asymptotic form of © 
when t - t Q -» oo is 


0 (k,p,t) « PKL^K[ e -vk 2 (t-to) . e -ttk 2 (t-t 0 )J (64) 

2(k - v) 2 


Finally, the quantity A' is of interest only in the combinat ion 
A + A 1 , which determines the Fourier transform of the correlation ■Suj/ 
(see eqs. (ll) and (12)). A consideration of the nondimens ional quantity 
(A + A 1 ) /Am shows, with the aid of equation (57) > that A 1 /Am = 

ofl /\lt - tA Hence, when t - t 0 


A + A 1 « A 


(65) 


where A is given by equation ( 65 ) . 

During the final period of decay, therefore, the asymptotic solutions 
do not depend on the initial conditions for the temperature field but only 
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on those for the velocity field, through the parameters Co and t Q . 
In the following sections, various properties of the temperature field 
will be derived on the basis of the approximate relations ( 63 ) to ( 65 ). 


Some Integral Properties of the Temperature Field 

There are several quantities, corresponding to mean values at one 
point only, that can be expressed entirely in terms of integrals involving 
the functi ons 0 and A. Before these results for the temperature 
field are presented, the corresponding results for the velocity field 
in terms of the energy spectrum function E will be listed for later use. 

From equation (40) and the definitions of reference 14, (pp. 47 
to 51 ) > 14 i 8 seen that the mean-square kinetic energy 1 b given by 


— , „ p°° 

u 2 = ^ uTu^ = - / E(k) dk = 

3 3 Jo 


_ \/at 


32 


v(t - to) 


- 5/2 


(66) 


the microscale A, by 


_ 1 _ 

A 2 


2 

15 u 2 





( 67 ) 


and the integral scale L, by 


L = 


2U 2 




( 68 ) 


Thus, in all further results it will be possible to express the quanti- 
ties C Q and t - to in terms of ( u 2 ) ^ and A, which are of more 
direct physical significance. The resulting formulas will then contain 
only quantities that are directly observable at any instant during the 
final period of decay, in addition to the given constant parameters of 
the problem (i.e., v, k, and (3). 
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Perhaps the quantity of greatest physical interest predicted by the 
present theory is the level of mean- square temperature fluctuations, which 
follows from equations (23), (64), (66), and (67) in the forms 



3\/2v2(i _ a ) 2 




2<y 


U + a 


1/2 


+ aV2 


1 g2 

5 (1 - a) 2 


h* 2 p 2 a 2 




+ a 1 / 2 


(69) 


The first form indicates that S 2 decays more 

slowly than vP- , just as Corrsin found in the simpler problem with constant 
mean temperature (ref. 9). 

Another quantity of some interest involves the mean-square nradient 
of d. Thus, if 


where 




1/2 


A/v. 


°i - r- (70) 

oXji * 

and if the components of G parallel ( Gp = A^G^ j and perpendicular to ~h 
are considered and axisymmetry noted, there results 

G 2 = G p 2 + 2G n 2 = f f k 2 0(k,p) dp dk (71) 

-1 

9 2 

The quantities Gp^ and G n are, respectively, the mean-square values 
of the component of. G parallel to A and the component in any direction 

A/ 

perpendicular to A and are given by 



NACA. TN hlQ6 


25 


G p 2 = k^9(t,l0 dp dk 


^L^C 0 [v(t-t 0 )]- 3/2 

60\/2v 2 (l - a} 2 




,3/2 


+ a3/2 


_ 1 _ g 2 
15 (1 - g) : 


R A 2p2 


1 - 2/'-22_? /2 + erS/ 2 


ll + ff 


( 72 ) 


and 




M? 


-) 8 (k,^) 


dp dk = 



( 73 ) 


From these results one can, of course, obtain the mean-square gradient 
of -8 in an arbitrary direction with respect to A . Relation ( 75 ) 

expresses the physical fact that the mean square of the molecular heat 
transfer fluctuations is greater "across the mean gradient" than it is 
"down the mean gradient," in a ratio that is independent of p . This 
conclusion is rather surprising, in view of the fact that the mean heat 
transfer (both molecular and turbulent) is entirely in the direction of 
the mean temperature gradient. 

Consider now the cross correlation L^(r,£). It is found, for 
example , that 


du p 


A ± Li(0,A) = at 






■)A(k,p) 


dp dk 


( 7*0 
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and 

dun = 0 (75) 

corresponding to the fact that the heat-transfer vector Li(0,A) is 
different from zero only in the direction of A. From equations (63) 
and (74) follows the explicit relation 



The heat-transfer correlation coefficient 



introduced by Corrsin (ref. l) now becomes^ 



^For isotropic turbulence, Up^ = u^. Also, note that |3 is by 
definition a positive constant. 
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A quantity of* particular physical significance is the ratio of 
mean turbulent heat transfer to mean molecular heat transfer given by 


-Sup 



2a 


il + ° 



(79) 


The behavior of this ratio is evidently just what would be expected phys- 
ically. For a given turbulent motion (i.e. , a given value of R 7 ), the 

ratio decreases with decreasing values of a (i.e., increasing thermal 
conductivity of the fluid) . On the other hand, for given fluid properties 
(i.e . , a given value of a), it increases with increasing Reynolds num- 
ber R^ of the turbulent motion. 


For the correlation Gjpj, where cdj = e , some rather 
Involved calculations lead to the relatively simple result 





e ijk Ak 


( 80 ) 


which is different from zero if and only if G and cu are perpendicular 
to each other as well as to ]\. After the general tensor relation for 

in terms of A is obtained, the above relation is derived most 

readily by ta k i n g rectangular coordinates x 2 , and with 

along A. It can then be shown that Gj_o>j_ = GgOVj = = 0 and, with 

the further use of spherical coordinates (k, 7 , and p) as described 

by the defining relations (4l) , that 


G ~ ~ 


= J' k2^A dk 


4ji 

3 



k^A dk 


(81) 


From these results and equation ( 63 ), equation (80) follows. 

Another quantity of some physical interest is the microscale or 
"dissipation" length parameter for the temperature field, first 
introduced by Corrsin (ref. 9) and defined by 



28 


NA.CA TN 4186 



»v 


In the present case* equation (82) reduces to 



(83) 


It should be noted that the geometrical interpretation of A,g is now 
not so simple as it was in the isotropic case (ref. 9)* In fact, with 
use of the relation m = r k A^ jr , it is found that 


d 2 0(r,A) _ d 2 0(r,m) 2 80(r,m) + 1 - m 2 de(r,m) 

8 r i ^ r 2 r dr dm r 2 dm 


(84) 


from which follows 



However, averaging this equation with respect to m leads to the simple 
result 



Thus, in an average sense, A$ is still a representative length scale 
for the variation of 9(r,m) in the vicinity of r = 0. 
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Another characteristic length for the temperature field is the 
thermal integral scale 


L* = 



(& 7 ) 


also introduced by Corrsin (ref. 9 ) • This quantity is expected to give 
a measure of the distance over which the temperature fluctuations are 
effectively correlated, as the integral scale L does for the velocity 
fluctuations. When 6 is expressed as the inverse of the transform rela- 
tion (22) and the spherical coordinates (k, 7 , and |3) in the k space 

described previously are used, equation ( 87 ) takes the form 


dp dk dr 

( 88 ) 


m l 

9(k,p) cos(krmp) J 
J 


kr\/(l - m 2 )(l - p 2 ) 


The integration with respect to r can be immediately carried out (see 
ref. 15 , p- 43), leading to 



When the expression ( 64 ) for 0 is introduced, the integration with 
respect to p is straightforward. The final integration with respect 
to k is readily performed by integrating from ko to » and then 
taking the limit as ko -* 0 , with the aid of the formula 
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log e x Q + 


poo 

2 / 
J x 0 \/A 


_ x 2 

xe logg x dx - 


2 log e 



dx 


(90) 


where x = k\j2v(t - t D ) . 
equation ( 69 ) , is 


The end result, after elimination of 


d 2 by 



+ a 


1/2 


-1 


log e 


(l + g) 2 
ka 


(91) 


An examination of the preceding results indicates that in this 
problem, for a given field of turbulence , many of the thermal character- 
istics are strongly dependent on the Prandtl number. It is instructive, 
therefore, to present a summary of the formulas derived above in the 
special cases: a » 1 (corresponding to several liquids, such as water 

and lubricating oils), a = 1 (approximately representative of gases), 
and a « 1 (corresponding to certain liquid metals). 

For a » 1, then, the following approximate relations result: 


=r - r x ( u 2) 1 px, 

- 2i . 1 R x 2 0 , 

3 A 


Gp 2 ar i .5/ VS 2 


Eh S - g ^--- x 0-1A 


& 

bii - 1 „ p(u?) 

G i“j = - — H >, - — r 

h * + 4 l08e 8 

1 -16 \ ' 


, 1/2 > ( 92 ) 

*X 


ql/2 



NACA TN 4186 


31 


Perhaps the most noteworthy feature here is that A^ and are "both 

much smaller than A, indicating that the temperature correlation curve 
is much more strongly curved near the origin than the velocity correlation 
curve and falls off to zero much more rapidly. It should he emphasized, of 
course, that this case places rather severe limitations on R^. Since the 

present theory is restricted to small values of P^ = (dA^/A^R^ = 


the Reynolds number R^ has to he extremely small if both the conditions 


cr » 1 and 


\j6dR- 


« 1 are to he satisfied. 


For a = 1, the formulas become 


'Y 



M _ 6 

A 5' 





In this case, the characteristic lengths for the temperature field are of 
the same order of magnitude as those for the velocity field. The result 
for Rh is very interesting, since it implies that the fluctuations -8 
and -Up = -Uj_A-j_ at a point are perfectly correlated. (Note that -Up 
is the velocity component in the direction of decreasing mean temperature.) 
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For. a « 1, the following approximate formulas are obtained: 


•a 2 S’ i o^R^p 2 ! 2 , 

3 


°p 2 - o 2 ^ 2 ? 2 


1/2 


au r 


S’ - I <%(u 2 ) pA, 


-au p _ i „ 2 o 
3f - R \ <S d . 

icp 3 A 


T> ~ ± 

h VS 


~ 1 


G i“j = - f ffR A 


0U 2 ) 


(94) 


6 ijk^k 


T- s6 ’ 


^,2f( ltn 2 )l08e i 


Although A^ is of the same order of magnitude as A in this case, 
is much larger, indicating that the temperature correlation function falls 
off rather slowly for large separation distances. In fact, the decrease 
of 9(r,m) for large values of r is so slow in the limiting case cr -» 0 
that the integral in equation ( 87 ) becomes divergent and becomes 

infinite. Thus, in this limiting situation the quantity loses its 

significance as a representative length for the temperature field. The 
length A-g remains significant, however, as a measure of the radius of 

curvature of the temperature correlation function near the origin. In 
the next section, a more detailed analysis of the behavior of the cor- 
relation functions for the temperature field will be presented. 

A final observation of considerable interest is that, if the arbitrary 
assumption is made initially that the turbulence is statistically station- 
ary^, then the final results agree with relations (94) when a is small. 
The above results for A^ and R^ were obtained by Mdlfese (ref. 8), 
and most of the others were obtained in a preliminary version of the 


^This implies that any mean value involving fluctuating quantities 
at different times is dependent only on the time differences and not on 
the individual times. Thus, mean values involving quantities at the same 
time, such as those in the present investigation, are independent of the 
time. The term "nondecaying" is also used in this sense. 
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present work, all on the basis of this hypothesis. Although stationary 
isotropic turbulence is impossible, strictly speaking, its assumption 
evidently leads to a formulation of the problem that is a valid approxi- 
mation when the Prandtl number is small. 

A clearer picture of the situation can be obtained as follows. 

First, it is noted that the assumption of stationary turbulence implies 
that the temperature field is also stationary, so that the time-derivative 
terms in equations ( 13 ) to ( 15 ) are absent and Li' = 0. On the other 
hand, if the assumption is not made, it is fo und that dQ/dt = 

0(Q/ (t - t 0 )) for any correlation Q, involving the temperature, so 
that in each of equations ( 13 ) and ( l4) the time-derivative term is 
0 [oKq j}<r) relative to the heat -conduct ion term and is t hus negligible 
when a is small. Moreover, it is also found that Li + Li' Li as 
t - to — > 00 • Hence, as far as the solution of the problem dur ing the 
final period of decay is concerned, the basic equations obtained with 
the assumption of s t at ionar it y are valid approximations for small values 
of cr, providing the correlation function f(r/A) describing the tur- 
bulence is properly specified. 

Physically, when the Prandtl number is small, the effects of time 
variation are negligible, and the dissipation through conduction is 
balanced primarily by the convective heat transfer associated with the 
mean temperature gradient. The dominating influence of the mean tem- 
perature gradient in this situation appears to be the reason for the 
basic difference between relation (94) for A^/a and the result obtained 
by Corrsin (ref. 9 ) to. the problem with zero mean temperature gradient, 

namely, that over a wide range of conditions A^/A = O^cr -1 / 2 ) for an 

Prandtl numbers. In the latter problem, the time variation is always 
a dominant factor. 


Some Detailed Properties of the Temperature Field 

The results of the preceding section were concerned with mean values 
at a single point and yielded only limited information about the various 
correlations as functions of the spatial separation. The characteristics 
of the temperature field associated with the two-point mean values are 
now considered in more detail. 

The scalar Fourier transform relation between 8(r,m) and ©(k,p) 
which follows from equation (22) is 
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nl 

®(k, (j.) 

j 


cos (knap) J Q 




dp dk 


(93) 


If a correlation coefficient for the temperature field is now defined in 
the usual way as 


M(r,m) = 

$ 2 


9 ( r , m ) 

9(0,m) 


(96) 


then using equations (64) and (95) gives 


M(y,m) 


IX * - 

•X 2 / 

2 _ e -x 2 /2ff 

1 eoa(xytim) J 0 

xy\J(l - m 2 )(l - [i 2 ) 

341 dx 

J 

X ( i • “ 2 > 

(e" x2 / 2 - e- x2 / 20 ) 2 dti ax 


(97) 


where y =» r j\j2v(t - t 0 ) = )[2r /'A. The integrals occurring in this 

expression can be evaluated quite easily (e.g., see ref. 15, p. 57, 
formula (52)) and the result written in the form 


M(y,m) = Mp(y) + m%2(y) 


(98) 


where 


Mi(y) = 


Mn(y) - 2) 


2a 


1+0 


M nm 


2 a 


'1 + a. 


\/oMn (y^) 


(99) 


2o 
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M2(y) = 


M 2 i(y) - 2 \J^~^ M 2i( y \Jj^~^) + ^ 2 i(y^) 


1 - 2 ' 


2a 


1 + a 


+ \[ a 


( 100 ) 


and 


Mn(y) = 


b y 


erf(y/2) + i 

y 


erf (y/2) 
y/2 


2 e -y2/4 

fit 


( 101 ) 


M2l(y) = r ^ i erf (y/2) - ^ 


4 y 


erf (y/2) _ _2_ e ~y 2 / 4 

y/2 


-g2 

e s d6. 


(102) 


The symbol erf(^) denotes the error function — 

fit J 0 

For small values of y, M(y,m) has the typical parabolic behavior 

^3/2 


l.JjS-) + a3/2 • 

„ (y , m) . x - i(i - f) — x ; g 1/2 ** + + ■*) 


i . 2^L_y v + cV 2 


(105) 


while for large values of y it decreases in the mann er 


-y 2 A _ 2 f 2g j 1//2 e -ffy 2 /2(l-Hi) + ff -l/2 e -ffy2/4 


M(y,m) 


3(1 - m2) 


1-2 




1 + 0 


1 + o r 
ay 2 


(104) 
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It Is noteworthy that, for 0 ^ 0, M(y,m) always becomes negative for 
sufficiently large values of y. 

The limiting form of M(y,m) for zero Prandtl number is of particular 
interest, in view of the remarks made in the last section. For 0 -» 0, in 
fact, 


M(y,m) = Mn(y) + m% 21 (y) 



(105) 


For small values of y, the limit function Mo(y,m) = Mn (y) + m^Mpi (y) 
has the parabolic behavior 


MoCy^m) 



(106) 


but for large values of y it decreases remarkably slowly, in the manner 




(107) 


The difference between the limit of relation (104) as a -> 0 and rela- 
tion (107) is due to the fact that the ratio (Mo - M)/Mq does not 
remain uniformly small for all values of y. For a very small but finite 
value of 0, the function M(y,m) follows the limit function Mo(y,m) 
very closely as y increases, but it begins to deviate appreciably when 
y becomes very large. In fact, It can be shown from equations (104) and 
(107) that Mo - M is of the order of magnitude of Mo when y is of 

the order of a' 1 / 2 . For » 1, |m| decreases to zero much more 

rapidly than Mq- The decrease of J M [ is, in fact, rapid enough to 

insure the convergence of the Integral of M (see eq. (87))* On the 
other hand, equation (107) shows that the integral of the limit func- 
tion Mo diverges logarithmically. These remarks help to explain the 
behavior of the thermal integral scale described previously. 

The case 0=1 is of considerable physical importance, since it 
corresponds approximately to gases, most of which have Prandtl numbers 
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slightly less than unity. The result for M(y,m) in this case, ■which 
is easily obtained by taking the limit a -> 1 in equations ( 98 ) to (100), 
has the especially simple form 


M(y,m) = e-y 2 A(l - 


The functions Mj_(y) and M 2 (y) have been evaluated numerically 

for the cases or = 0 , 0 . 01 , 0 . 1 , 0 - 71 , 1 , 10 , and 100 , respectively. 

The results for all other cases can be readily obtained from those for 
a = 0 in table i4, since the functions M]_(y) and Mp(y) for this case 
are the functions M i 1 (y) and Mpi fy) occurring in the general for- 
mulas (99) to (102). The numerical values for the case a = O.'Jl, cor- 
responding to air under standard laboratory conditions, are presented in 
table 3 , since they may be of seme general interest. 

The resultant correlations M(y,m) for m = 0 and 1 (corresponding 
to a = 90 ° and 0 °, respectively) are shown in figures l(a) to l(c), with 

the velocity correlation f(y) = e“^A included in each case for com- 
parison. The correlations for all other values of m lie between those 
shown. The curves for a = 0.01 and 0.1 are included in figure 1 along 
with those for a = 0 in order to illustrate the approach to the limit 
discussed previously. In contrast to the slow rate of decrease of M(y,m) 
for small values of a, the rapid rate of decrease for large values of a 
is to be noted. This behavior, it may be recalled, was indicated by the 
previous results for the thermal integral scale. The physical significance 
is just that for small values of the Prandtl number, the temperature fluc- 
tuations remain significantly correlated over much greater distances than 
the velocity fluctuations, while for large values the opposite is true. 

Consider now the correlation L^(r,A). It Is observed that if 


■6Up' - Aj_Lj_^r, A'j (3-09) 


d) + y2 e -y 2 A 
4/ 4 


( 108 ) 


A'here is a loss of accuracy near a = 1 when eqs. (99) and ( 100 ) 
are used directly, with functional values interpolated from table 1 . 
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(see fig. 2(a)), then the Fourier transform (eq. (24)) yields 


*Up‘ - i 1 ~ ^ 2 )k 2 A(k,ti) cos(krm|i) 

(no) 

Thus, letting 


■Sup' 

N(r,m) = — — (111) 

i5up 


kr 


\j{ 1 - m 2 )(i - n 2 ) 


d|a dk 


and substituting for A(k,p.) from equation (63) give 


a ( x - 


N(y,m) 


1+tJ. 


-x2 


- e 2a i coa > 




(xypm) J Q xy\j(x - m 2 )(l - p 2 ) 


dp. dx 


(13J2) 


dp dx 


where x 
to give 


:\j2v(t - to) and y = r /\/ 2v (t - to). 


This can be evaluated 


N(y,m) = Nq(y) + m2N 2 (y) 


(113) 


where 


%(y) 


H n(y) - 


I 2 a 
[l + a 


5/2 


N 


1 - 


2a 

,1 + 


nly' 

3/2 


2a 

1 + a 


(114) 
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and 



z< 

e-y 2 /*» - £ 

erf (y/2) 2 r _ v 2 fh 

> 

2 

1 3,2 

y/2 fit 



N 2 i(y) = - I 


s -y2 jk _ 3Vjc 


erf(y/2) _ _2_ e _ y 2 jk 

y/2 |/5t 


For small values of y, N(y,m) also has the familiar parabolic 

> 5/2 


N(y, 


■,*) - 1 - ifl - 2^LAL2 lSL y2 + oUl + a*) 

10 V 2/ , _ ^3/2 ' ' 




3/2 


while for large values of y it decreases in the manner 



Just as in the case of M(y,m), for cr ^ 0, the function N(y,m) 
becomes negative for sufficiently large values of y. 


(115) 


(116) 


(HT) 

behavior 


(118) 


(119) 

always 
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The limiting form of N(y,m) for small Prandtl numbers is also of 
interest. For a -* 0 the following relation is obtained: 

N(y,m) = N^Cy) + m 2 ^!^) + o((j5/ 2 ) (120) 

For small values of y, the limit function N 0 (y,m) = Nn(y) + m^N p-| (y) 
has the behavior 

N 0 (y,m) -• 1 - - ^jy 2 + 0(y4) . (121) 

For large values of y and for m 4 l/^3 it decreases somewhat faster" 
than M 0 (y,m) , in the manner 

— +o(e-y 2 

y 5 V 


The case - m = l/y3, for which 


N 0 (y>l/\/3) = e“ y2 / 4 (123) 

divides the functions N 0 (y,m) into those that are positive everywhere 
and those that become negative (see fig. 3(a)). 5 The differences between 
the behavior of the limit function N 0 (y,m) for large values of y and 
that of the actual function N(y,m), for any small but finite value of cr, 
are to noted. The reasons for these differences are very much like those 
described for the function M(y,m) . 




5 it might be noted that there is no such dividing correlation when 
a 4 0, since then the function N(y,m) always becomes negative for 
large values of y. 
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In contrast to M(y,m), the correlation N(y,m) 
for large Prandtl numbers. Thus, -when a -> 


has a limiting form 


N(y,m) 


(2^~- lj _1 2\f2N 11 (y^2) - N^y) 



(124) 


The limit function in this case is obviously not radically different in 
shape and asymptotic behavior from the function N(y,m) for a value of a 
of order unity. 

For ff a 1, a simple result is again obtained, in fact, exactly the 
same functional form for N(y,m) as for M(y,m), namely. 


N(y,m) 



(125) 


The functions Np(y) and N 2 (y) have been computed for the cases 

a = 0, 0.1, 0.71, 1, 10, and «. The numerical values for the cases of 
most interest, that is, a = 0 and a = 0 . 71 , are tabulated in tables 2 
and 3, respectively. The resultant correlations for m = 0 and 1 (cor- 
responding to a = 90° and 0°) are shown in figures 3( a ) to 3(c) . The 

velocity correlation f(y) = e -y2 /^ is also included in each figure. 

For a = 0, f(y) is the same as (see eq.. ( 123 ) ) . 

From these results for the correlation flup' it is possible to derive 
the corresponding results for the defining scalars of L±^r,A^ and hence 

to obtain the correlation of ■& with the velocity in an arbitrary direc- 
tion. Thus, from equations (l7), (19), and (20), 


\iI*i^r,A) = (l - n^rDr + 2 L(r,m) = 3up Nj_(y) + m% 2 (y) 


( 126 ) 


an equation that can only be satisfied provided that L(r,m) is inde- 
pendent of m, in which case 
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L(y) = | 3u p [Ni(y) + N2(y) 


(127) 


From these results, it follows that the correlation ■3u n ' (see fig. 2(b)) 
is given by 


— — = m^l - m2Ng(y) =» — sin 2a Ng(y) (128) 

3u p 2 


The function N2(y) is plotted in figure 4 for cr = 0, 0.1, 0.J1, 1, 10, 
and oo. 

Since the correlation of •S with the component of u-j_ orthogonal 
to the plane formed by r and A vanishes by axisymmetry, it is necessary 
to consider only the correlation of -S with a velocity vector coplanar 
with r and A. The results given here for dUp’ and tfu^ then suf- 
fice to determine directly the correlation of S with the velocity in 
an arbitrary direction without recourse to its scalar representation, 
such as those shown, for exanple, in figures 2(c) and 2(d). 

The analytical and numerical results for N(y,m) described above 
Indicate that, although N decreases rather slowly for large values 
of y when a is small, the Prandtl number dependence is not other- 
wise very strong. This is in marked contrast to the results for M(y,m), 
especially those for large Prandtl numbers . In particular, although 
M(y,m) decreases at a much more rapid rate than the velocity correla- 
tion f(y) when the Prandtl number is large, N(y,m) still decreases 
at a rate not significantly different from that for f(y), just as it 
does for Prandtl numbers of order unity. 

Mathematically, these conclusions are to some extent suggested by 
the forms of the basic equations (15) and (l4) for the double temperature 
correlation 0 and the temperature -velocity correlation L±, respectively. 

In particular, equation (13) involves only the parameter k, while equa- 
tion (14) involves both k and v in the combination k + v. The form 

of equation (l4) then indicates that, when the Prandtl number is large 
and k + v ~ v, Lj_ is approximately independent of the Prandtl number 
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and has a characteristic length scale essentially the same as that for 
R-^j (i.e., the microscale A). On the other hand, for small Prandtl 

numbers, k + v s k, and Li might then be expected to depend on the 

Prandtl number. Finally, equation (13) always involves k explicitly 
and k + v implicitly (in the term involving Li), so that 0 might be 

expected to depend on k and thus on the Prandtl number in general. 

Physically, the temperature -velocity correlation is governed by both 
momentum transfer and heat transfer, with only the convective part of the 
latter (associated with the mean temperature gradient) being appreciable 
when the thermal conductivity is small. Thus, it is not too surprising 
that for small thermal conductivity this correlation behaves much like 
the velocity correlation. The double temperature correlation, on the 
other hand, depends only on heat transfer, and in this case the molecular 
part is always as important as the convective part. Consequently, the 
conclusions for the temperature correlation are also to some extent 
expected. Consider, for example, two fluids in turbulent motion with 
essentially the same dynamical and physical properties, except that one 
has low thermal conductivity (large Prandtl number) and the other high 
thermal conductivity (small Prandtl number) . In order for the tempera- 
ture fluctuations at two points to follow each other at all closely (and 
thus be significantly correlated), it would be expected that the two points 
would have to be rather close together in the first case but that they 
could be much farther apart in the second. It should be recalled, of 
course, that the present investigation is restricted to low Reynolds 
numbers, when molecular transport processes are the controlling factors. 

For high Reynolds numbers, the inertial transport processes neglected 
here are usually more important, and the final results may be quite dif- 
ferent from the present ones. 


Temperature and Heat-Transfer Spectrum Functions 

In the last two sections the characteristics of the temperature field 
revealed by the one-point and two-point mean values have been studied. 

From some viewpoints, a further gain in physical insight is provided by 
an examination of the properties of the Fourier transforms of the correla- 
tion functions . These transforms and some functions closely related to 
them will now be considered in more detail. 

Since the behavior of the general three-dimensional Fourier transform 
may in general be rather complicated, it is often convenient to define a 
simpler function by averaging over all directions of the vector argument k. 

Thus, by analogy with the definition of the energy- spectrum function for the 
velocity field, the three-dimensional (or spherically averaged) temperature- 
spectrum function may be defined as 
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0(k) = 



(129) 


This is just the integral of the Fourier transform ©(k, p.)/2rtk2 of the 

7 A 

correlation M over the surface of a sphere of radius k. Thus, 0 

is the density of contributions to $2 associated with wave numbers of 

magnitude k, so that 



( 130 ) 


From equations (64) and (129) the following results are obtained: 


©(k) = i\fl E A ^ 2 A 3 (e" x2 / 2 - e~ x2 f 2a ) (ljl) 

5 (1 - a) 2 


where x 


k^ 2v(t - t 0 ) . 


The maximum value of 0(k) 


is given by 



a 1 " 0 R a % 2 A 5 


(132) 


and occurs at 


x m ~ 



(133) 


The function ©(k) / ^ has been evaluated numerically for Prandtl 

numbers a = 0, 0.01, 0.1, 0.71, 1, 10-, 100, and », and the results are 
shown in figure 5* The numerical values for the case a = 0.71 are 
given in table 4. The wide range of variation of the results with the 
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Prandtl number is to be noted, as well as the do min ance of the small 
wave-number contributions (associated with the large-scale structure of 
the temperature field) for small Prandtl numbers and the dominance of 
the large wave-number contributions (associated with the small-scale 
structure) for large Prandtl numbers. This behavior of the spectrum 
function is perhaps to be expected, on the basis of the previous con- 
clusions for the correlation function. 

~ k 

The behavior 0 ~ if for small values of the wave number, which is 
the same as it is for the turbulence -energy spectrim E(k), is to be com- 
pared with Corrsin's result 0 « k.2 in the problem with constant mean 
temperature (ref. 16) . It should be kept in mind that under considera- 
tion is the asymptotic part of the spectrum for large values of the time, 
which has characteristics determined mainly by the mean temperature gra- 
dient and the turbulence . Even in the present problem a part of the 
spectrum beginning with k^ could be present, in general, as Mdlfese 
(ref. 8) has shown (see eqs . ( 38 ) and ( 55) ) • However, aB demonstrated 
previously, even if such a part were present initially, it would even- 
tually make only a negligible contribution to the total spectrum function. 

In a similar manner, a three-dimensional spectrum fun ction A(k) 
is now defined for the mean turbulent heat transfer -dup, where the 

negative sign is chosen in order to have a positive quantity (i.e., 

-•3up = | 'flu!, as a result of the definition of A) : 


A(k) = -2rck 2 



dp 


(134) 


Thus, A(k) is the density of contributions to -dup associated with 
wave numbers of magnitude k, and 




dk 


(135) 
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Equations (30), (65), and (134) then give the results 


A(k) = 


- ^ k 2 A(k) 


■ Ilf rh R x(“ 2 ) 1/2 ^ 2x2 ( v e ' : 




x 2 _ e 2a 


(136) 


The maximum value Am occurs at a value of x given by 


- 


„ 2 , 
x m 


1 - e 


2a 


1 -a 


1 - L-f-g e 2tJ 

2 a 


-xm - 1 


(137) 


The function A(k)/Am has been evaluated numerically for Prandtl 
numbers 0 = 0, 0.1, O.71, 1, 10, and «>, and the results are shown in 
figure 6. The numerical values for the case 0 = 0.7I are given in 

table 4. In contrast to the conclusions for ©(k)/^, the variation of 

A.(k)/Am with the Prandtl number is not large. In fact, the main con- 
tributions always tend to come from essentially the same range of wave 
numbers as that corresponding to the energy- containing eddies of the 
turbulence . 

Equation (136) shows that A(k) begins with a term of order k^ 
for small values of k during the final period of decay, just as it 
does, in general, initially (see eqs. (36) and (56))- The general con- 
clusion that the heat-transfer spectrum begins with at least the fourth 
power of k was obtained previously by M^lfese (ref. 8). 

Although the three-dimensional spectrum, as described in the pre- 
ceding two cases, has perhaps the greatest physical significance, it 
cannot be directly measured. One-dimensional spectrum functions are 
more suitable for experimental purposes, since they can be obtained 
fairly easily from correlation measurements and in fact can sometimes 
be measured directly. Some useful results relating such one-dimensional 
spectrum functions with the three-dimensional ones are now presented. 
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A one -dimens io nal Fourier transform of the double temperature cor- 
relation 0^r, A^ = - 6 - 6 ' can be made along a line at an angle aq with A 

as follows: 


9 I 



( 2 k )- 1 J ©(^,0,0) e - ^ 11 * 1 drq 


(138) 


where the dependence on mq = cos ctq due to the axisymmetry of the tem- 
perature field has been explicitly indicated. Since 0^r, Aj = 0(r,m) is 
an even function of both r and m, equation ( 158 ) can be written as 





(139) 


Then, with the aid of equation (95 ) t this becomes 


= f SqJqS 0 008 ( k i r i) cos (^l 10 !^) J o - m i 2 ) - v- 2 ) ^ 


(H*o) 


After substitution for ©(k,p) in terms of ©(k) from equation (131), 
the integration with respect to p and rq can be easily carried out 

in the cases mq = 0 and mq = 1 (see ref. 15 , p. 45), with the results 


e l( k l»°) 




(141) 


( 142 ) 
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Since 0 m-^ = 0^r-]_,O^ + m^ 2 9^r-j_,l^ - 

is seen to be 


, the general formula 


0 l(kl' m l) = (l - mi 2 ) 0 l( k l^o) + m 1 a Sl I ( k 1 A) 


(145) 


The inverse relations to these are as follows: 


§(k) = - i k 
5 




'(k,o) - r 

^0 


k 0i'( k ,o) 


dk 


(144) 


_d_ 

©l'(k,l)' 

dk 

k 


e(k) = - 


4k 


3(l - 


2m-\ 2 


m l^ + l 


(145) 


01 '^m-^ - k 1 " 111 ! 2 J k m l dk 


(146) 


The prime denotes differentiation with respect to k. The integral from 
0 to k in relations (l44) and (l46) can be replaced by the integral 
from oo to k, since it can be shown that the integral from 0 to oo 
is zero in each case. From any one of these relations, then, the three- 
dimensional temperature spectrum can be . determined from a measured one- 
dimensional spectrum. It appears that formula (l44) should be the most 
useful since 0i(hi>o) is a Fourier transform in a direction normal to A. 

In an experimental study of the present problem, the mean flow direction 
would likely be normal to A, so that a direct measurement of 0j^k]_,O) 

would be possible (subject to certain approximations), without having to 
determine it from a correlation measurement. 

It might be mentioned that in the limit a -* 0 the one -dimens ional 
spectrum function ©j is well behaved for small values of k, in spite 

of an apparent logarithmic divergence of^the integrals in equations (l4l) 
and (142) in the limiting case. Since 0 has a maximum value of the 

order of p2 a t a value of x of the order of (a log e a)^"/ 2 (see 
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equations ( 1 J 2 ) and (133)), it follows from equations (l4l) and (l42) 
that 0j has a maximum value of the order of cr^logg cr. Although this 

A 

is much larger than the maximum value of 0 , it still approaches zero in 
the limit cr -»■ 0 . 

For the temperature-velocity correlation Lj_(r,A) = since it is 

a vector quantity, many one -dimensional spectrum functions could be defined. 
Two of the most useful ones appear to be the following: 


Ai^kp^mJ = (2n)' 1 J Lp(r 1 ,0,0je _lk l r l dr-L 

i r 

= * J 0 L p( r i jm i) cos ( k i r i) 

A I i(ki,m 1 ) = Ln^O^Oje-^l ^ 

= “ f Q COB ( k l r l) ^1 


(147) 


(148) 


With the aid of equations (110) and ( 136 ) , the following special results 
can be obtained in the manner described previously: 


Ai(k 17 0) 



(149) 


Ac (^1*1) 



( 150 ) 


From these, since Lp^r-^m^ = Lp^r-j_,0j + m^ |%> ( r ^, 1 j - Lp^yOj 
general case follows: 


the 
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Ai(ki,m 1 ) = (l - m 1 2)A I (k 1 ,0) + m^'k^l) 


(15D 


Also, since Ln^r^m^ = m-^1 - m-j_ 2 Lp^r^l) - Lp^r^ojj , (see fig. 2(b)), 
the result for Ajj is 


A ii(^i* m i) = ^ly 1 - m i 2 A i( k iA) - A i( k i>°) 


(152) 


The inverse relations to relations (l49), (150), and (152) are 


A(k) = | k 

5 


A I 


•(k,o) - r 

d 0 


k A l’(k,0) 


dk 


(155) 


A(k) = - — k5 — 
3 dk 


A l'(k,l) 


( 154 ) 


A(k) = - 


4k 


^ly 1 " m i L 


A II * (k, n^i) - ^ kAjx’^^m^ dk 


(155) 


The integrals from 0 to k may again be replaced by integrals from °° 
to k. 


Although the one -dimensional spectrum Aj^k]_,oj could perhaps be' 

measured directly in the experimental situation, such a determination is 
not so easy as that for 9x[k-j_,0), and a more satisfactory procedure 


appears to be to obtain a one-dimensional spectrum from a measured cor- 
relation, such as L n ^r^,m 1 ^, for example. From the corresponding one- 


dimensional spectrum A n(kx,mx), the three-dimensional heat-transfer 
spectrum A.(k) could then be obtained by means of equation ( 155) • 
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It should, be emphasized that all of the special relations between 
the various spectrum functions derived above are valid only under the 
special restrictions of the present investigation. 


CONCLUDING REMARKS 


In the present investigation of heat transfer in isotropic turbu- 
lence with a constant mean temperature gradient a variety of special 
results have been derived for the final period of decay. From these 
results the following major conclusions emerge: 

(a) During the final period of decay, the temperature field becomes 
asymptotically independent of the initial conditions on the temperature. 

Its characteristics eventually depend essenti ally only on the mean tem- 
perature gradient, the physical properties of the fluid, and the char- 
acteristics of the turbulence. 

(b) For a given field of turbulence, a.1 1 results depending on tem- 
perature fluctuations alone are strongly dependent on the Prandtl number. 

Examples are the results for the mean square -S 2 of the temperature 
fluctuations and for the double-temperature-correlation coefficient 

57* 'M It is particularly noteworthy that for small Prandtl numbers the 
temperature fluctuations remain significantly correlated over much larger 
distances than the velocity fluctuations, while for large Prandtl numbers 
the opposite is true. 

(c) Many of the results depending on both the temperature and velocity 
fluctuations display much less dependence on the Prandtl number. For 

example, the mean turbulence heat transfer •Sup is of the same order of 

magnitude for very large Prandtl numbers as it is for a Prandtl number 
of unity. The Prandtl number variation of such quantitie s is usually 
more important for small Prandtl numbers. The quantity 6Up, for example, 

approaches zero when the Prandtl number approaches zero. However, the 
behavior of the temperature -velocity correlation coefficient 3up ' j -Sup 

changes relatively little over the full range of Prandtl numbers, in 
marked contrast to the behavior of the double temperature correlation 
coefficient . 

(d) Corresponding trends are present in the behavior of the three- 
dimensional spectrum functions. For the temperature spectrum, the sm all 
wave-number part is most important when the Prandtl number is small, while 
the large wave-number part is most Imp ort ant, when the Prandtl number is 
large. For the heat-transfer spectrum, however, the main contributions 
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always come from essentially the same range of wave numbers as that 
corresponding to the energy-containing eddies of the turbulence. 

In evaluating these conclusions, the restriction to small Reynolds 
and Pdclet numbers should be kept in mind; however, although the results 
are strictly valid only in the limit of vanishing R^ and any 

radical changes for small but finite values of these parameters would 
not be expected. Experience with other such_approximations suggests 
that the present results might even be quantitatively accurate for values 
of these parameters up to 1 and perhaps qualitatively correct up to values 
of about 5* It would, of course, be desirable to investigate the devia- 
tion from the limiting results for zero R^ and by means of a 

second-order approximation, in which terms of order R^ and are 

retained. In such an analysis, the triple correlations, which are neg- 
lected in the present study, would have to be considered. Thus, in 
addition to corrections to the present results, conclusions regarding 
the role of the triple correlations for small Reynolds and Pdclet numbers, 
which are of great interest in themselves, would also be obtained. 

The situation at high Reynolds and Pdclet numbers, it should be 
emphasized, is likely to be quite different from that described by the 
present theory or any extension of it in the manner just outlined. In 
these circumstances, the inertial transfer terms in the basic equations, 
which should be small for the present type of analysis to be valid, are 
normally dominant. Consequently, a different mathematical approach has 
to be used, and the results would be expected to be qualitatively quite 
different. For example, one would not expect to find a strong Prandtl 
number dependence of the properties associated with the large-scale 
structure of the temperature and velocity fields, as found here in some 
cases, although such a dependence would likely still be present in the 
characteristics of the small-scale structure. In an investigation of 
this situation, the statistical hypothesis mentioned in the introduction 
appears to be essential in formulating a mathematical problem that is at 
all tractable. As suggested previously, there are still serious analytical 
difficulties, but it may be expected that eventually this approach will 
lead to a solution of the problem. 


The Johns Hopkins University, 

Baltimore, Md. , February 25, 1957> 

and 

Ballistic Research Laboratories, U. S. Army, 

Aberdeen Proving Ground, Md., February 25, 1957* 
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TABLE 1.- FUNCTIONS DESCRIBING TEMPERATURE 
CORRELATION M(y,m) = Mx(y) + m^Cy) 
FOR FRANDTL NUMBER a = 0 


y 

Mi(y) 

M2(y) 

Mi(y) + M2(y) 

0 

1.0000 

0 

1.0000 

.2 

.9960 

.0020 

.9980 

.4 

.9842 

.0079 

• 9921 

.6 

.9650 

.0173 

.9823 

.8 

• 9392 

.02 99 

.9691 

1.0 

.9076 

.0450 

.9526 


.8713 

.0619 

• 9332 

■ 

.8315 

.0800 

• 9114 


.7893 

.0984 

.8877 

1.8 

.7459 

.1165 

.8624 

2.0 

.7022 

.1338 

.8360 


.6592 

.1498 

.8090 


.6176 

.1641 

.7817 


• 5779 

.1765 

.7544 

2.8 

.5405 

.1870 

.7275 

3.0 

.5057 

.1954 

.7011 

3.2 

.4735 

.2019 

.6754 

3-4 

.4440 

.20 66 

.6306 

3.6 

.4171 

.2097 

.6268 

3-8 

.3926 

.2114 

.6040 

4.0 

• 3704 

.2119 

• 5823 

4.2 

.3503 

.2113 

.5616 

4.4 

• 3321 

.2099 

.5420 

4.6 

.3156 

.2079 

.5234 

4.8 

.3006 

.2053 

.5059 

5-0 

.2869 

.2023 

.4892 

5-2 

.2744 

.1991 

, .4735 

5.4 

.2630 

.1956 

.4586 

5.6 

.2525 

.1920 

.4445 

5.8 

.2428 

.1883 

.4311 

6.0 

.2339 

.1846 

.4185 

6.5 

.2142 

.1755 

.3897 

7-0 

•1977 

.1667 

.3643 

7-5 

.1835 

.1583 

.3419 

8.0 

.1714 

.1306 

■ 3219 

8.5 

.1607 

.1434 

.3041 

9-0 

.1514 

.1368 

.2881 

9- 5 

.1430 

.1306 

• 2737 

10.0 

.1356 

.1250 

.2606 
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TABLE 2.- FUNCTIONS DESCRIBING TEMFERATURE-VELOCITI 
CORRELATION N(y,m) = N x (y) + m2N 2 (y) 

FOR FRANDTL NUMBER a - 0 


y 

Ni(y) 

N 2 (y) 

Ni(y) + N2(y) 

0 

1.0000 

0 

1.0000 

.2 

.9881 

.0058 


.4 

.9550 

.0233 


.6 

.8970 

.0507 


.8 

.8256 

.0857 

• 9093 

1.0 

.7569 

.1256 

. 8626 

1.2 

.6418 

.1675 

.8094 

1.4 

.5452 

.2083 

.7515 

1.6 

.4455 

.2455 

.6909 

1.8 

.5526 

.2768 

.6294 

2.0 

.2676 

.3008 

.5684 

2.2 

.1926 

.3169 

.5094 

2.4 

.1286 

• 3249 

.4535 

2.6 

.0761 

.3253 

.4014 

2.8 

.0545 

• 3190 

• 3535 

5-0 

.0030 

.3073 

.3103 

5-2 

-.0198 

.2914 

.2716 

5-4 

-.0555 

.2727 

.2373 

5.6 

-.0449 

.2522 

.2071 

5.8 

-.0500 

.2312 

.1812 

4.0 

-.0518 

.2103 

.1585 

4.2 

-.0515 

.1902 

• 1390 

4.4 

-.0492 

.1714 

.1222 

4.6 

-.0463 

.1540 

.1077 

4.8 

-.0429 

.1382 

.0953 

5.0 

-.0394 

.1240 

.0846 

5.2 

-.0559 

.1113 

.0754 

5.4 

-.0327 

.1001 

.0674 

5.6 

-.02 96 

.0901 

.0605 

5-8 

-.0269 

.0814 

.0545 

6.0 

-.0244 

.0736 

.0492 

6.5 

-.0195 

.0580 

.0387 

7.0 

-.0155 

.0465 

.0310 

* 7-5 

-.0126 

• 0378 

.0252 

8.0 

-.0104 

.0312 

.0208 

8.5 

-.0087 

.0260 

•0173 

9.0 

-.0075 

.0219 

.0146 

9-5 

-.0062 

.0186 

.0125 

10.0 

-.0053 

.0160 

.0106 
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TABLE 3.- FUNCTIONS DESCRIBING TEMPERATURE CORRELATION 
M(y,m) = Mi(y) + m^fe(y) AND TEMEERATURE-VELOCITy 
CORRELATION N(y) = Ni(y) + m%2(y) FOR 
ERANDTL NUMBER a = O. 7 I 


y 

Ml 

M2 

Mi + M2 

*L 

»2 

N x + N 2 

0 

1.000 

0 

1.000 

1.000 

0 

1.000 

.4 

.954 

.033 

.967 

.929 

•035 

.964 

.8 

•755 

.115 

.870 

.737 

.126 

.863 

1.2 

• 514 

.218 

• 752 

.482 

• 237 

•719 

1.6 

.270 

.309 

•579 

.231 

.326 

•557 

2.0 

.071 

• 359 

.430 

.034 

.366 

.400 

2.4 

-.060 

• 358 

.298 

-.085 

•352 

.267 

2.8 

-.124 

.317 

• 193 

-.131 

.298 

.167 

5-2 

-.155 

.248 

.115 

-.128 

.225 

.097 

3.6 

-.110 

.178 

.068 

-.101 

•155 

.052 

4.0 

-.084 

.120 

.036 

-.069 

•095 

.026 

5-0 

-.024 

.024 

.000 

-.016 

.019 

.003 
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TABLE 4 .- TEMPERATURE AND HEAT-TRANSEER SPECTRUM 
FUNCTIONS 8 (x)/ 0 m AND A(x)/Am FOR 
ERANDTL NUMBER cr = 0-71 


X 

A f, S 

@/©m 

/s /s 

A/^Cl 

0 

0 

0 

.2 

.0041 

.0034 

.4 

.0561 

.0479 

.6 

.2229 

.1944 

.8 

.5035 

.4516 

1.0 

.7983 

.7422 

1.2 

.9784 

.9496 

1.4 

• 9751 

.9947 

1.6 

.8151 

.8801 

1.8 

.5837 

.6709 

2.0 

.3631 

.4468 

2.2 

.1983 

.2625 

2.4 

.0959 

.1571 

2.8 

.0160 

.0269 

5-2 

.0017 

.0035 

3.6 

.0001 

.0003 

4.0 

.0000 

.0000 
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Figure 1.- Concluded 


m = I . a - 0° 

m-0 , a ~ 90° 

f (y) 



KACA TO 4186 
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Qin I 

£■ = N L (y) + m% 2 (y) • 

0Up 



(t) = | sin 2a N 2 (y) . 

0up 2 



= cos a 




(4) 



sin a N-j_(y) . 


Figure 2.- Temperature-velocity correlations. 



0 



(a) Prandtl numbers a = 0 and cr =■» 0.10. 

Figure 3.- Temperature-velocity correlation coefficient 

N(y,m) s Qn'^/eu^. m = cos otj y = I/2" r/Xj X a Turbulence microscale^ 

jr f r vi 

f (y) = Velocity correlation. 


991-1] iCL VOW. 




z 


3 


(c) Prandtl numbe 
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Figure 4.- Temperature-velocity correlation coefficient 
N 2 (y) = (2/Bin 2a)eu , n /e^ for various Prandtl numbers . 
f(y) = Velocity correlation; y ■» l[2 r/?v; = Turbulence microscale. 
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Figure 5 .- Kuree-dimensional temperature-spectrmn function 9(x) far 
various Prandtl numbers, x ~ kA/'/2; 6^ = Maximum value Of 8) 

Tv =* Turbulence microscale: E(x) = Turbulence energy spectrum. 
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Figure 6.- Three-dimensional heat-transfer spectrum function A(x) for 
various Prandtl numbers . x = lcX/^2; ^ « Maximum value of 

X =. Turbulence microscale; E(x) = Turbulence energy spectrum. 
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